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Abstract.  A  South  China  Sea  warm  pool  with  sea  surface  temperature  (SST) 
higher  than  29.5°C,  recently  reported  by  Chu  and  Chang  [1995a,  b]  and  Chu  et  al. 
[1997],  appears  in  the  central  South  China  Sea  (west  of  the  Luzon  Island)  in  boreal 
spring,  strengthens  until  the  onset  of  the  summer  monsoon  (mid-May),  and  then 
weakens  and  disappears  at  the  end  of  May.  The  transient  features  and  interannual 
variabilities  of  the  warm  pool  have  not  yet  been  studied.  Here  we  use  a  subset  of  the 
U.S.  Navy’s  Master  Oceanographic  Observation  Data  Set  (MOODS)  to  investigate 
the  surface  thermal  features.  First,  we  employed  an  optimal  interpolation  scheme  to 
build  up  a  10-day  interval  synoptic  data  set  for  December  1963  to  November  1984  on 
a  0.5°  x  1°  grids  (finer  resolution  in  zonal  direction)  from  the  MOODS  SST  data.  An 
ensemble  mean  SST  field  (T)  was  established  with  a  rather  weak  horizontal  gradient 
(28.5°C  near  the  Palawan  Island  to  26°C  near  the  southeast  China  coast).  Second, 
we  performed  a  composite  analysis  to  obtain  the  averaged  SST  anomaly  field  T 
deviating  from  the  ensemble  mean  for  the  winter  and  spring  seasons  (December- 
May).  During  December-March,  T  is  negative  almost  everywhere  throughout  the 
whole  South  China  Sea.  In  early  April,  positive  T  with  closed  isoline  (warm  pool) 
was  evident  west  of  Luzon  Island.  In  May,  the  central  SCS  warm  anomaly  becomes 
stronger.  On  May  11-20,  the  central  SCS  warm  pool  (1140-119°E,  14°-19°N) 
has  T  >  1.8°C.  The  size  of  the  warm  pool  is  around  200,000  km2.  Third,  we 
performed  an  empirical  orthogonal  function  (EOF)  analysis  on  the  residue  data 

(T),  deviating  from  T  +  T,  for  the  winter  and  spring  seasons,  in  order  to  obtain 
transient  and  interannual  variations  of  the  SST  fields.  EOF1  accounts  for  35.5%  of 
the  variance  and  resembles  the  ensemble  mean  pattern  of  nearly  parallel  contours 
with  a  maximum  value  in  the  southeast  and  a  minimum  value  in  the  northwest. 
EOF2  accounts  for  21.4%  of  the  variance  and  is  characterized  by  a  warm/cool  pool 
(116°-118°E,  16°-18°N)  west  of  the  Luzon  Island.  The  corresponding  principal 
component  {PC2)  has  strong  interannual  variability  with  a  maximum  value  of  10  on 
February  11-20,  1965  and  a  minimum  of  -12  on  March  21-31,  1964.  This  indicates 
the  appearance  of  either  a  warm  pool  with  a  maximum  strength  of  1.2°C  or  a  cool 
pool  with  a  maximum  strength  of -1°C.  Combination  of  f  and  PC2x EOF2  leads  to 
an  occurrence  of  a  central  SCS  warm  pool  from  April  to  May  with  a  warm  anomaly 
varying  between  0.8°  and  3°C. 


1.  Introduction 

The  South  China  Sea  (SCS)  has  a  bottom  topography 
(Figure  1)  that  makes  it  a  unique  semi-enclosed  ocean 
basin  that  is  overlaid  by  a  pronounced  monsoon  sur- 
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face  wind.  Extended  continental  shelves  (less  than  100 
m  deep)  are  found  on  the  western  and  southern  parts, 
while  steep  slopes  with  almost  no  shelves  are  found  in 
the  eastern  part  of  SCS.  The  deepest  water  is  confined 
to  a  bowl-type  trench.  The  maximum  depth  is  around 
4700  m. 

On  the  basis  of  limited  data  sets,  both  cool  and  warm 
pools  exist  in  SCS.  Dale  [1959]  reported  a  cool  pool  off 
the  central  Vietnamese  coast  in  summer.  Nitain  [1970] 
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Figure  1.  Geography  and  isobaths  showing  the  bottom 
topography  of  the  South  China  Sea. 

found  a  cool  pool  located  to  the  northwest  of  Luzon. 
Reports  from  the  South  China  Sea  Institute  of  Oceanol¬ 
ogy  (SCSIO)  [SCSIO,  1985]  indicate  that  in  the  central 
South  China  Sea,  a  warm  pool  appears  in  both  sum¬ 
mer  and  winter  but  closer  to  Vietnam  in  summer  at 
the  surface.  Recently,  a  warm  pool  was  reported  in 
the  central  SCS  during  the  late  spring  season  [Chu  and 
Chang,  1995a, b;  Chu  et  al.,  1997;  Tseng,  1995]  and  a 
cool  pool  was  detected  in  the  central  SCS  during  De¬ 
cember  29,  1993  to  January  5,  1994,  from  the  analysis 
of  TOPEX/POSEIDON  data  [Soong,  et  ah,  1995], 

An  international  South  China  Sea  Monsoon  Exper¬ 
iment  (SCSMEX)  was  recently  initiated.  Two  of  the 
major  objectives  are  to  identify  the  influence  of  heat¬ 
ing  contrasts  between  the  SCS  and  surrounding  regions 
and  the  role  of  early  monsoon  (April-May)  convection 
and  multiscale  processes  in  the  SCS  in  the  abrupt  tran¬ 
sition  and  subsequent  evolution  of  the  East  Asian  mon¬ 
soon  [SCSMEX  Science  Plan,  1995].  As  a  part  of  the 
SCSMEX  effort,  our  interest  was  concentrated  on  the 
spatial  and  temporal  variabilities  of  the  SCS  thermal 
fields  from  December  to  May. 

2.  Master  Oceanographic  Observational 
Data  Set  (MOODS) 

We  use  the  extensive  U.S.  Navy’s  Master  Observa¬ 
tional  Oceanographic  Data  Set  (MOODS)  for  the  study. 
The  MOODS  is  a  compilation  of  ocean  data  observed 
worldwide  consisting  of  (1)  temperature-only  profiles, 


(2)  both  temperature  and  salinity  profiles,  (3)  sound- 
speed  profiles,  and  (4)  surface  temperature  (drifting 
buoy).  These  measurements  are,  in  general,  irregular 
in  time  and  space.  Due  to  the  shear  size  (more  than 
six  million  profiles  total  for  the  global  ocean)  and  con¬ 
stant  influx  of  data  to  the  Naval  Oceanographic  Office 
(NAVOCEANO)  from  various  sources,  quality  control 
is  very  important.  The  primary  editing  procedure  in¬ 
cluded  removal  of  profiles  with  obviously  erroneous  lo¬ 
cation,  profiles  with  large  spikes  (temperature  higher 
than  35°C  and  lower  than  -2°do  not  match  the  charac¬ 
teristics  of  surrounding  profiles,  such  as  profiles  showing 
increase  of  temperature  with  depth.  After  quality  con¬ 
trol,  the  historical  MOODS  data  (1900-1995)  consisted 
of  189,059  profiles  for  the  whole  SCS  (5°S-25°N,  105°- 
120°  E).  Our  study  domain  is  a  subarea  of  SCS  covering 
9°  to  25°N  and  109°to  120° E,  and  our  investigation  pe¬ 
riod  is  December  1963  to  November  1984.  The  total 
number  of  profiles  of  this  subset  is  70,546.  For  con¬ 
venience,  we  group  the  previous  year’s  December  data 
with  the  current  year’s  January-November  data;  for  ex¬ 
ample,  the  1964  data  set  includes  the  data  of  December 
1963  to  November  1964. 

The  temporal  and  spatial  distribution  of  MOODS 
data  is  irregular.  Certain  periods  and  areas  are  very 
well  sampled,  while  others  lack  enough  observations  to 
gain  any  meaningful  insights.  There  exist  some  10  to  20- 
day  gaps  with  no  observations  in  the  whole  SCS.  Figure 
2  shows  the  sparsity  of  profiles  in  the  southeast  portion 
of  the  study  domain  (from  Liyue  Bank,  Zhenghe  Reef 
to  Yingqing  Reef)  and  in  the  coastal  region  of  China 
continent.  Figure  3  indicates  a  heavily  sampled  period 
during  the  Vietnamese  War  (1965-1969).  The  maxi¬ 
mum  number  of  observations  is  in  1965  (around  7200 
profiles).  The  minimum  number  of  observations  is  in 
1984  (near  980  profiles).  Furthermore,  vertical  resolu¬ 
tion  and  data  quality  are  also  highly  variable  depending 
much  on  instrument  type  and  sampling  expertise.  Tem¬ 
poral  and  spatial  irregularities  along  with  the  data  res¬ 
olution  and  quality  problems  must  be  carefully  weighed 
in  order  to  avoid  mathematically  induced  variability. 

3.  Establishment  of  Gridded  Data 

We  use  1°  x  1°  monthly  sea  surface  temperature  (SST) 
climatology  [Levitus,  1982]  as  the  mean  field  and  an 
optimum  interpolation  scheme  to  obtain  gridded  data 
from  the  MOODS  data.  For  1964-1984,  the  SST  per¬ 
turbations  for  SST  observations  were  ob¬ 

tained  by  subtracting  the  climatological  value  at  the 
closest  grid  from  the  observation.  Optimum  interpola¬ 
tion  assigns  a  weight  to  each  observation  that  accounts 
for  variation  in  spatial  sampling.  The  interpolated  tem¬ 
perature  anomaly  at  the  grid  point,  T'G ,  is  a  linear  com¬ 
bination  of  the  observed  anomalies,  T[,  T], ...,  T'N  with 
weights  ai,  «2,  •••)  <xn'- 

T'g  (1) 

i= 1 


CHU  ET  AL.:  SOUTH  CHINA  SEA  WARM  POOL  DETECTED  IN  SPRING 


15,763 


MOODS  Stations  for  64-84  (Total:70546) 


r 


4 


0 

*D 


3 


110  115  120 

Longitude  (E) 

Figure  2.  Spatial  distribution  of  MOODS  stations  during  1964-1984  for  this  study. 


The  values  of  the  weights,  a,(t  =  1, 2, TV),  are  found 
by  minimizing,  in  a  least  square  sense,  the  difference 
between  the  interpolated  value  at  the  grid  point  and 
the  true  value  there  and  are  obtained  from  solving  the 
algebraic  equations  [Gandin,  1965]: 

N 

^2  +  A_2a,  =  HGi,  *  =  1,2, ...,  TV.  (2) 

3  =  1 

Here  A  is  the  signal-to-noise  ratio,  mj  is  the  autocorre¬ 
lation  between  locations  i  and  j ,  and  nGi  is  the  auto¬ 
correlation  between  the  grid  point  and  location  i.  The 
autocorrelation  function  (ACF)  is  defined  by 

1  fR 

A*  =  ^2  y  T1  (r0)T' (r0  +  r)dr0  (3) 

where  r0  denotes  the  independent  space  vector  defining 
the  location  of  points  in  a  sampling  space  R,  r  is  the 
space  lag,  and  s2  is  the  variance.  Here  fi  is  computed 
by  paring  the  anomalies  into  bins  depending  upon  their 
separation  in  space,  r.  The  values  of  p  will  be  obtained 
from  calculating  the  correlation  coefficient  for  all  the 


anomaly  pairs  in  each  bin  and  represented  by  a  Gaus¬ 
sian  model  for  the  combination  of  different  lags 

1 l(m)  =  r?(0)  exp  [— A2(mAr)2]  (4) 

where  r/(m)  denotes  the  ACF  value  in  the  bin  with  the 
horizontal  separation  mAr.  Ar  =10  km  is  the  incre¬ 
ment  for  the  space/time  separation.  is  the  horizon¬ 
tal  decorrelation  scale.  Therefore  the  autocorrelations 
l-i i j  and  (IQ,  are  given  by 

Hij  =  r?(m0),  p.Gi  =  T}{mGi )  (5) 

where  m,j  and  mG,  represent  spatial  separations  be¬ 
tween  locations  i  and  j  and  between  location  i  and  grid 
point,  respectively. 

We  applied  this  optimal  interpolation  method  to  build 
up  a  near- 10-day  interval  synoptic  data  set  from  Decem¬ 
ber  1963  to  November  1984  on  a  0.5°  x  1°  grid  (finer 
resolution  in  the  zonal  direction)  from  the  MOODS  tem¬ 
perature  data.  Our  total  gridded  data  set  consists  of 
756  consecutive  SST  fields. 
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Figure  3.  Temporal  distribution  of  MOODS  stations  during  1964-1984  for  this  study. 


4.  Composite  Analysis 


T(xi,yj,it)  =T(xi,yj,ti)-T{xi,yj)  (8) 


We  now  examine  the  data  to  see  if  we  can  obtain 
a  warm  pool  signal  in  the  mean  seasonal  {winter  to 
spring)  data.  After  the  gridded  data  set  is  established, 
the  SST  data  are  represented  by  T(xi,  t/j,  Tfe,  i(),  where 
(xi,yj)  is  the  horizontal  grids,  t*  =  1964,  1965,  ..., 
1984  is  the  time  sequence  in  years,  and  tj  =1,  2,  ...,  36 
is  the  time  sequence  within  a  year  in  10-day  intervals. 
The  first  10-day  interval  in  the  year  starts  from  Decem¬ 
ber;  that  is,  ti  —  1  represents  December  1-10  (previous 
year),  1/  =  2  denotes  December  11-20  (previous  year), 
...,  ti  =  4  refers  to  January  1-10  (current  year),  and 
ti  =  36  indicates  November  21-30  (current  year).  Be¬ 
fore  investigating  the  annual  variation  of  SST,  we  define 
the  following  two  temporal  averages: 

T(xj, yj,ti)  =  ^ ^  'y^T[xj,yj,Tk,ti) 
k 

At  =  21(1984  —  1964)  (6) 

which  are  the  annual  mean  values,  and 


36 

1  V — '  =, 


T(xi,  j/j )  _  y  ^  T[xj,  yj,  ti) 


which  are  the  ensemble  mean  values.  The  ensemble 
mean  (1964-1984)  SST  field  over  the  SCS  (Figure  4) 
shows  a  pattern  of  northeast-southwest  oriented  isotherms 
with  a  positive  temperature  gradient  toward  southeast 
near  the  equator.  The  ensemble  mean  has  a  rather 
weak  horizontal  temperature  gradient.  T  decreases 
from  28.5°C  west  of  Palawan  to  26°C  near  the  southeast 
China  coast. 

The  annual  mean  values  deviated  from  the  ensemble 
mean,  T(x,,  j/j), 


lead  to  the  composite  features  of  the  10-day  mean 
SST  anomalies.  Figures  5  and  6  show  T(x,-,  t/j, fj), 
t;  =  1, 2, ...,  18,  since  we  are  only  interested  in  the  ther¬ 
mal  fields  for  the  first  part  of  the  year  (December  to 
May), 

The  10-day  mean  SST  anomaly  (T)  has  the  following 
features: 

1.  The  first  feature  is  a  cooling  phase.  During 
December  to  March,  T  is  negative  almost  everywhere 
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Figure  4.  The  ensemble  mean  of  the  SST  field  during 
1964-1984. 
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Figure  5.  Averaged  10-day  SST  anomalies  relative  to  the  ensemble  mean  from  December  to 
February:  (a)  December  1-10,  (b)  December  11-20,  (c)  December  21-31,  (d)  January  1-10  fe) 
January  11-20,  (f)  January  21-31,  (g)  February  1-10,  (h)  February  11-20,  (i)  February  21-28  (or 
29).  Solid  (dashed)  contours  denote  positive  (negative)  anomalies.  v 


throughout  the  whole  domain  and  becomes  more  neg¬ 
ative  from  the  west  coast  of  Palawan  to  the  southeast 
coast  of  China.  This  indicates  that  the  SCS  cooling 
during  the  winter  monsoon  season  is  the  weakest  in  the 
southeast  SCS  and  the  strongest  near  the  Taiwan  Strait 
(<  — 3°C). 

2.  The  second  feature  is  an  early  warming  phase.  In 
early  April,  positive  T  is  evident  west  of  Luzon  Island 
and  west  of  Palawan  Island.  At  the  Luzon  Strait,  T 


stays  negative  until  late  April  (April  21-30).  This  means 
that  early  SCS  warming,  starting  in  the  southeast  SCS 
and  advecting  to  the  central  SCS,  is  not  forced  by  the 
Kurushio  intrusion.  It  also  implies  that  the  Kuroshio 
water  entering  SCS  in  April  might  be  cooler  than  the 
SCS  water  (at  least  at  the  surface). 

3.  The  third  stage  is  a  mature  warming  phase.  In 
May,  the  central  SCS  warm  pool  becomes  stronger  than 
the  early  warming  phase.  During  May  11-20,  the  central 


15,766 


CHU  ET  AL.:  SOUTH  CHINA  SEA  WARM  POOL  DETECTED  IN  SPRING 


Figure  6.  Averaged  10-day  SST  anomalies  relative  to  the  ensemble  mean  from  March  to  May: 
(a)  March  1-10,  (b)  March  11-20,  (c)  March  21-31,  (d)  April  1-10,  (e)  April  11-20,  (f)  April  21-30, 
(g)  May  1-10,  (h)  May  11-20,  (i)  May  21-31.  Solid  (dashed)  contours  denote  positive  (negative) 
anomalies. 


SCS  warm  pool  (114°-119°E,  14°-19°N)  has  f  >  1.8°C 
and  accounts  for  the  maximum  mean  warm  anomaly. 
The  size  of  the  warm  pool  is  around  200,000  km2. 

5.  Empirical  Orthogonal  Function 
(EOF)  Analysis 

The  composite  analysis  shows  the  existence  of  a  cen¬ 
tral  SCS  warm  pool  in  the  10-day  mean  SST  anomaly 
data  during  April  to  May.  What  are  the  transient  fea¬ 
tures  of  the  SCS  warm  pool?  What  are  the  synoptic 
and  interannual  variabilities  of  SCS  warm  pool?  We 


used  empirical  orthogonal  function  (EOF)  analysis  to 
investigate  these  issues. 

5.1.  EOFs  for  the  SST  Anomalies 
The  SST  synoptic  anomalies  obtained  by 

T{xi,yjTkti )  =  T{xi,yj,Tktl)  -  T{xityjtt )  (9) 

are  re-arranged  into  a  N  x  P  matrix,  T(rn,tp),n  = 
1,2,...,  A;  and  p  =  1,2 ,  ...,P.  Here  P  =  378  is  the 
total  number  of  time  points  used  for  computing  the  co- 
variance  matrix,  i.e.,  21  years  of  every  10-day  measure- 
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ments  from  December  to  May;  N  =  255  corresponds 
to  the  number  of  grids  ( i  —  1,2,  ..,17;  j  =  1,2,, ..15). 
Empirical  orthogonal  function  analysis  [Lorenz,  1956], 
widely  used  in  oceanographic  and  meteorological  re¬ 
search  [e.g.,  Weare  et  al.,  1976]  (see  review  by  Rich- 
man,  [1986])  is  the  same  as  principal  component  (PC) 
analysis  [Hotelling,  1933]  in  the  statistics  community. 
PCs  are  the  amplitudes,  which  are  functions  of  time,  of 
their  corresponding  EOFs.  These  EOFs  can  be  found 
by  calculating  the  unitary  eigenvectors  of  the  covariance 
matrix  associated  with  the  sample  data  field.  EOF  anal¬ 
ysis  separates  the  data  sets  into  eigenmodes.  Generally 
speaking,  each  mode  has  an  associated  variance,  nondi- 
mensional  spatial  pattern,  and  dimensional  time  series. 
One  nice  feature  of  the  EOF  analysis  is  that  it  does  not 
require  the  data  to  be  continuous  in  time.  We  can  ana¬ 
lyze  the  data  from  December  to  May  to  investigate  the 
variability  during  that  period.  Therefore  we  have  378 
time  points  in  the  EOF  analysis.  From  this  data  matrix 
a  255-square  spatial  covariance  matrix  is  calculated  by 


■ftu 

Hi  2 

R\n 

R  = 

R21 

R22 

R2N 

Rn  1 

RN2  - 

-  Rnn 

Rnm  —  y/ p  )  ' T(r„, <p)T(rm , tp), 

N  —  255,  P  =  378  (10) 

where  n  and  m  (1,2,  ...,N)  denote  the  grid  locations. 
The  diagonal  elements  of  the  covariance  matrix  Rnn  (n  = 

1,2,  ...,1V)  are  the  variance  at  location  rn.  The  off- 
diagonal  elements  are  the  covariance  with  spatial  lag 
equal  to  the  difference  between  the  row  and  column  in¬ 
dices.  This  symmetric  matrix  has  its  N  real  eigenvalues 
Xa  and  eigenvectors  <f>a(r j),  such  that 

N 

^  —  Ani^Q,(r,),  i—  1,2,  ...,1V  (11) 

1=1 

The  eigenvectors  <f>i,4>2t—,<f>N  are  called  empirical  or¬ 
thogonal  functions.  Each  <j>a  is  a  255-point  (17x15  grid 
in  this  study)  distribution  of  SST  anomaly  pattern.  The 
eigenvalues,  Aq(q  =  1,2,  ...,1V),  are  all  positive,  and  the 
summation  of  them,  A a ,  equals  the  total  variance. 
Therefore  \a  is  considered  as  the  portion  of  total  vari¬ 
ance  “explained”  by  the  EOF  <f>a.  It  is  convenient  to 
label  the  eigenfunctions  4>a  so  that  the  eigenvalues  are 
in  descending  order,  i.e., 

Ai>A2>A3>...  (12) 

5.2.  Principal  EOF  Modes 

Intraseasonal  and  interannual  anomalies  can  be  quan¬ 
titatively  investigated  with  the  method  of  EOF  analysis 
and  a  number  of  its  generalized  forms.  In  order  to  delin¬ 
eate  the  major  modes  of  variability  in  winter  and  spring 


seasons  (December  to  May)  of  21  years  (1964-1984),  we 
perform  an  EOF  analysis  using  the  method  outlined  in 
section  5.1  and  obtain  the  first  six  leading  EOFs,  which 
are  able  to  account  for  almost  81.7%  of  the  total  vari¬ 
ance  during  the  seasons  (Table  1).  When  combined, 
EOF1  and  EOF2  explain  over  half  (57%)  of  the  total 
variance.  Each  EOF  mode  is  normalized  so  that  its  to¬ 
tal  spatial  variance  is  equal  to  unity.  So  those  patterns 
of  the  first  six  EOFs  are  enough  to  explain  the  spatial 
anomalies  of  whole  SCS.  Hence  it  suffices  to  focus  on  the 
first  six  EOFs.  Their  structures  are  shown  in  Figure  7. 

The  EOF1  mode  (Figure  7a)  has  a  NE-SW  orienta¬ 
tion  and  accounts  for  up  to  35.6%  of  the  total  spatial 
variance  during  December  to  May.  This  pattern  is  be¬ 
lieved  to  be  related  to  the  ensemble  mean  pattern.  This 
indicates  that  such  a  pattern  (NE-SW  orientation)  re¬ 
peats  in  time.  To  illustrate  the  possibility  of  the  same 
pattern  occurring  in  the  mean  and  EOF1  fields,  we  as¬ 
sume  that  a  data  set  contains  only  one  pattern  A  with 
the  amplitude  varying  in  time.  Then  the  mean  field  is 
simply  the  pattern  A  with  the  average  amplitude,  and 
the  EOF  (only  one)  is  also  the  pattern  A.  By  removing 
the  mean  from  the  data  field,  the  average  value  of  the 
EOF  amplitudes  now  vanishes.  The  EOF2  mode  (Fig¬ 
ure  7b)  has  a  closed  isoline  of -0.1  to  the  west  of  Luzon 
Island  (116a-118°E,  16°-18°N)  and  accounts  for  about 
21.4%  of  variance.  This  warm/cool  pool  is  embedded 
into  the  warm  pool  identified  by  the  composite  analysis. 
However,  we  still  cannot  identify  whether  EOF2  repre¬ 
sents  a  warm  or  cool  pool  since  the  values  in  Figure  7 
are  normalized.  We  should  consider  the  product  of  the 
EOF  values  and  the  corresponding  time  series  of  the 
amplitudes,  also  called  PC  components. 

5.3.  Temporal  Variabilities 

The  data  matrix,  T(xi,yj,Tk,ti ),  is  thus  written  by 

Tixuyj^Ji)  =  Y.PC«&p)M*i,yj)  (13) 

a 

where  PCa(tp)  is  the  principal  component  with  a  unit  of 
degree  Celsius  and  a  size  of  P,  representing  the  tempo¬ 
ral  variation  of  the  associated  spatial  pattern  described 
by  EOF  <j>a(xi,yj).  The  time  series  is  analogous  to  a 
projection  of  SST  anomaly  through  the  “filter”  of  an 
EOF  mode  during  a  time  scale. 


Table  1. 

Variances  of  the  First  Six  Leading  EOFs 

EOF 

Variance,  % 

Cumulative  Variance,  % 

1 

35.6 

35.6 

2 

21.4 

57.0 

3 

12.4 

69.4 

4 

5.3 

74.7 

5 

4.7 

79.4 

6 

2.3 

81.7 
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Figure  7.  First  six  EOF  modes  (the  values  were  multiplied  by  10)  computed  from  January-May 
data  during  1964-1984:  (a)  EOF1,  (b)  EOF2,  (c)  EOF3,  (d)  EOF4,  (e)  EOF5,  and  (f)  EOF6. 
Solid  (dashed)  contours  denote  positive  (negative)  values. 


The  first  principal  component,  F‘Ci(fp),  for  1964- 
1984  is  shown  in  Figure  8.  The  EOF1  mode  cj>i(xi,yj) 
is  always  negative  throughout  the  whole  domain  (Fig¬ 
ure  7a).  Therefore  PCi(tp)  >  0  corresponds  to  nega¬ 
tive  SST  anomalies,  and  PCi(tp)  <  0  corresponds  to 
positive  SST  anomalies.  Taking  the  year  of  1968  as 
an  example,  PCi(tp)  is  always  positive  from  early  De¬ 
cember  (1967)  to  late  April  (1968),  which  implies  the 
negative  SST  anomalies  in  the  whole  SCS  during  this 
period.  Furthermore,  the  variability  of  EOF1  (Figure 
7a)  is  small  west  of  Palawan  Island  and  increases  north¬ 
westward  toward  the  southeast  coast  of  China. 

The  second  principal  component,  PC-2(tp)>  for  1964- 
1984  is  shown  in  Figure  9.  The  EOF2  mode  <j>2{xi,yj) 
is  negative  throughout  the  whole  domain  and  shows  a 


closed  isoline  of  -0.1  to  the  west  of  Luzon  Island  (Fig¬ 
ure  7b).  As  before,  PCi(tp)  >  0  corresponds  to  neg¬ 
ative  SST  anomalies,  and  PCiiip)  <  0  corresponds  to 
positive  SST  anomalies.  Taking  the  year  of  1964  as 
an  example,  PC2(tp)  is  almost  always  negative  from 
December  1963  to  early  May  and  reaches  a  minimum 
value  of  -12  on  March  21-31,  which  implies  positive  SST 
anomalies  over  the  whole  SCS  with  a  warm  pool  com¬ 
ponent  to  the  west  of  Luzon  Island  during  1964  winter 
and  spring  seasons.  The  maximum  SST  anomaly  of  the 
warm  pool  in  1964  is  around  1.2°C  [—0.1  x  (—12°)].  On 
the  other  hand,  taking  the  year  of  1965  as  another  ex¬ 
ample,  PC2(tp)  has  a  maximum  value  of  10  (mid  Febru¬ 
ary),  which  implies  negative  SST  anomalies  in  the  whole 
domain  and  a  cool  pool  component  to  the  west  of  Lu- 
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Figure  8.  Time  series  of  PCX  (December  to  May)  for  1964-1984.  The  unit  is  degrees  Celsius. 

zon  Island  with  a  magnitude  of  -1°C  (-0.1xl0°C)  on  warm  pool  formation  in  spring.  When  the  surface  wind 
February  11-20,  1965.  This  indicates  the  appearance  of  stress  curl  over  the  central  SCS  is  anticyclonic,  Ekman 
either  a  warm  pool  with  a  maximum  strength  of  1.2°C  downwelling  will  occur  in  the  central  part  of  the  bowl 
or  a  cool  pool  with  a  maximum  strength  of  -1°C.  and  upwelling  will  occur  near  the  boundary  of  the  bowl 

Combination  of  T  and  PC2XEOF2  leads  to  an  oc-  through  mass  balance.  The  downwelling  prevents  the 
currence  of  a  central  SCS  warm  pool  from  April  to  May  deep  cold  water  from  advecting  upward,  and  the  up- 
with  a  warm  anomaly  varying  between  0.8°  and  3°C.  welling  helps  the  deep  cold  water  advecting  upward. 

This  makes  the  central  part  of  the  bowl  warm  and  the 
6.  A  Possible  Combined  side  Part  °f  the  bowl  cool  and  leads  to  the  generation 

Wind- Topography  Effect  of  the  central  SCS  warm  pool. 

From  late  winter  to  early  spring,  a  surface  anticyclone 
The  bowl-type  bottom  topography  shown  in  Figure  usually  appears  over  the  central  SCS  [Cheang,  1987], 
1  provides  a  favorable  condition  for  the  central  SCS  This  anticyclone  generates  downwelling  in  the  central 
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Figure  9.  Time  series  of  PC2  (December  to  May)  for  1964-1984.  The  unit  is  degrees  Celsius. 


SCS  and  in  turn  prevents  the  cold  deep  water  from  be¬ 
ing  advected  to  the  surface.  As  spring  starts  the  cold 
northeast  (winter)  monsoon  diminishes  and  the  sky  over 
the  SCS  also  enters  a  more  clear  period  with  less  cloud 
cover  and  rainfall  and  rapidly  increasing  solar  radiative 
warming  at  the  sea  surface.  This  warm  pool,  with  a 
temperature  (between  0.8°  and  3°C)  higher  than  the 
surroundings  in  a  generally  warm  SCS,  may  lower  the 
atmospheric  surface  pressure  and  promote  the  onset  of 
the  southwest  (summer)  monsoon.  The  appearance  of 
a  cyclonic  circulation  with  the  lower  pressure  after  the 
summer  monsoon  onset  may  generate  Ekman  upwelling 
in  the  central  SCS,  bringing  the  deep  cold  water  into 


the  surface  mixed  layer.  This  process  may  then  destroy 
the  central  SCS  warm  pool. 

7.  Conclusions 

Our  study  shows  the  existence  of  a  central  SCS  warm 
pool  during  winter  and  spring  seasons  from  analysis  of 
the  Navy’s  Master  Observational  Oceanographic  Data 
Set  (MOODS).  Optimal  interpolation,  composite  anal¬ 
ysis,  and  EOF  analysis  were  used.  We  obtained  the 
following  results: 

1.  The  ensemble  mean  SST  field  (T)  was  established 
with  a  rather  weak  horizontal  gradient  (28.5°C  west 
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of  Palawan  Island  to  26°  C  near  the  southeast  China 
Coast). 

2.  The  composite  analysis  indicates  the  following  fea¬ 
tures  of  the  averaged  SST  anomaly  fields  T  relative  to 
the  ensemble  mean  in  the  winter  and  spring  ^seasons 
(December-May).  During  December  to  March,  T  is  neg¬ 
ative  almost  everywhere  throughout  the  whole  domain. 

In  early  April,  positive  T  with  closed  isoline  (warm 
pool)  is  found  west  of  Luzon  Island.  The  size  of  the 
pool  is  around  200,000  km2.  In  May,  the  central  SCS 
warm  anomaly  becomes  significant.  On  May  11-20,  the 
central  SCS  warm  pool  has  T  >  1.8°C. 

3.  EOF  analysis  was  applied  to  the  residual  SST 
data  T  (deviated  from  T  +  T)  to  investigate  the  tran¬ 
sient  and  interannual  SST  variabilities.  EOF1  accounts 
for  35.5%  of  the  variance  and  resembles  the  ensemble 
mean  pattern  of  nearly  parallel  contours  with  a  max¬ 
imum  value  in  the  southeast  and  a  minimum  value  in 
the  northwest.  EOF2  accounts  for  21.4%  of  the  vari¬ 
ance  and  is  characterized  by  a  warm/cool  pool  (em¬ 
bedded  into  the  warm  pool  identified  by  the  composite 
analysis)  to  west  of  Luzon  Island.  The  corresponding 
principal  component  ( PC2 )  has  strong  interannual  vari¬ 
ability  with  a  maximum  value  of  10  on  February  11-20, 
1965  and  a  minimum  of -12  on  March  21-31,  1964.  This 
leads  to  the  appearance  of  either  a  warm  pool  with  a 
maximum  strength  of  1.2°C  or  a  cool  pool  with  a  max¬ 
imum  strength  of-l°C.  In  other  words,  PC2XEOF2  is 
bounded  between  -1°  and  1.2°C. 

4.  Combination  of  T  and  PC2XEOF2  leads  to  an 
occurrence  of  a  central  SCS  warm  pool  from  April  to 
May  with  a  warm  anomaly  varying  between  0.8°  and 
3°C. 

5.  A  combined  wind-topography  effect  is  proposed 
to  explain  the  formation  and  destruction  of  the  central 
SCS  warm  pool.  On  the  other  hand,  the  appearance 
and  disappearance  of  the  central  SCS  warm  pool  will  af¬ 
fect  the  atmospheric  circulation.  Thus  the  central  SCS 
warm  pool  might  be  an  important  component  for  the 
SCS  coupled  air-ocean  system. 
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